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ABSTRACT 

We have measured the central structural properties for a sample of SO-Sbc galaxies down to scales 
of ~10 pc using Hubble Space Telescope NICMOS images. We find that the photometric masses of the 
central star clusters, which occur in 58% of our sample, are related to their host bulge masses such 
that Mps = 10 7 - 75±015 (Mbui/10 10 M q ) 076±0 - 13 . Put together with recent data on bulges hosting 
supermassive black holes, we infer a non-linear dependency of the 'Central Massive Object' mass on 
the host bulge mass such that M CM o/M Q = 10 7 - 51±0 06 (A4 B ui/10 10 X ) a84±0 06 . We argue that 
the linear relation presented by Ferrarese et al. is biased at the low-mass end by the inclusion of 
the disc light from lenticular galaxies in their sample. Matching our NICMOS data with wider-field, 
ground-based -fT-band images enabled us to sample from the nucleus to the disk-dominated region of 
each galaxy, and thus to perform a proper bulge-disk decomposition. We found that the majority of 
our galaxies (~90%) possess central light excesses which can be modeled with an inner exponential 
and/or an unresolved point source in the case of the nuclear star clusters. All the extended nuclear 
components, with sizes of a few hundred pc, have disky isophotes, which suggest that they may be 
inner disks, rings, or bars; their colors are redder than those of the underlying bulge, arguing against 
a recent origin for their stellar populations. Surface brightness profiles (of the total galaxy light, and 
the bulge component on its own) rise inward to the resolution limit of the data, with a continuous 
distribution of logarithmic slopes from the low values typical of dwarf ellipticals (0.1 < 7 < 0.3) to 
the high values (7 ~ 1) typical of intermediate luminosity ellipticals; the nuclear slope bi-modality 
reported by others is not present in our sample. 

Subject headings: galaxies : spiral — galaxies : structure — galaxies : nuclei 



1. INTRODUCTION 

The Hubble Space Telescope (HST) enables the study of 
the inner regions of nearby bulges and ellipticals down to 
spatial scales of ~10 pc, roughly one order of magnitude 
closer to the center than is feasible with typical ground- 
based data. These inner regions contain a small fraction 
of the ellipsoid mass, but they harbor the highest density 
regions of the galaxies and contain useful clues to their 
formation. 

The availability of NIR array detectors in the nineties 
fostered significant progress in the understanding of 
many aspects of bulges, including structural parameters, 
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colors, dust content and stellar populations, as well as the 
scaling of disk and bulge parameters, using ground-based 
imaging (e.g. Andredakis, Peletier, & Balcells 1995, here- 
after APB95; de Jong 1996; Seigar & James 1998; 
Knapen et al. 1995; Khosroshahi et al. 2000; Graham 
2001a, hereafter G01; Mollenhoff & Heidt 2001; Graham 
2001b, 2002; Eskridge et al. 2002; MacArthur, Courteau, 
& Holtzman 2003; Castro- Rodriguez & Garzon 2003). 

NIR data helped to establish that exponential pro- 
files provide better fits to the surface brightness pro- 
files of bulges than 14 models (Kent et al. 1991 for the 
MW bulge; Andredakis & Sanders 1994; de Jong 1996), 
and soon thereafter it was demonstrated that profiles of 
bulges of all Hubble types admit a particularly simple 
fit using the Sersic (1963; see Graham & Driver 2005) 
function 

I(R) = 1(0) exp{-6 n (R/R c ) 1/n } (1) 

(APB95; G01; Mollenhoff & Heidt 2001; MacArthur et 
al. 2003; see Caon et al. 2003 for the case of elliptical 
galaxies). In eqn. [1] R c is the half-light radius of the 
bulge, and b n w 1.9992n - 0.3271. The Sersic index n, 
which measures the curvature of the surface brightness 
profile, scales with bulge-to-disk luminosity ratio (B/D) 
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and with bulge luminosity. The Sersic index also provides 
a concentration parameter (Trujillo et al. 2001) which 
strongly correlates with the velocity dispersion and cen- 
tral supermassive black hole mass (Graham et al. 2001a, 
2001b), hence it is linked to global physical parameters 
of the spheroid. Numerical simulations also suggest that 
bulges have a range of profile shapes. Aguerri, Balcells, 
& Peletier (2001) show that the accretion of dense satel- 
lites onto disk-bulge-halo galaxies yields a growth of both 
the Sersic index and B/D, hinting that n may be linked 
to the accretion history and to the growth of bulges. A 
cold dark matter cosmological simulations of galaxy for- 
mation yield bulge-disk structures where the bulge pro- 
file shape ranges from exponential to 14 (Scannapieco & 
Tissera 2003; Sommer-Larsen, Goth, & Portinari 2003). 

The results given above, derived from ground-based 
data, bear the question of whether the inner regions to 
which the HST gives access also follow the Sersic func- 
tion. Our picture of elliptical galaxy nuclei had to be re- 
vised in several ways after the HST imaging campaigns. 
Giant ellipticals often show a rather sudden inward flat- 
tening of their surface brightness profiles, confirming the 
result from ground-based data that some ellipticals have 
"cores" (Kormendy 1985), while intermediate- luminosity 
ellipticals (—18 < Mb < —20.5) do not show cores; their 
profiles approach power laws throughout the inner re- 
gions ("power-law" galaxies); see Faber et al. (1997, here- 
after F97) and Rest et al. (2001, hereafter R01). Inner 
profile slopes decrease toward fainter luminosities, and, 
for dwarf ellipticals, approach the slopes seen in the nu- 
clei of giant, core galaxies, although dwarfs do not show 
profile discontinuities, i.e., do not show 'cores' (Graham 
& Guzman 2003; Ferrarese et al. 2006, hereafter F06). 
Many cores of ellipticals and SOs are dusty, and a frac- 
tion of them harbor central unresolved sources at HST 
resolution (Lauer et al. 1995, hereafter L95; Phillips et al. 
1996; Carollo et al. 1997; Ravindranath et al. 2001; Sti- 
avelli et al. 2001). Inasmuch as bulges share global sim- 
ilarities with ellipticals when studied from the ground, 
we enquire whether bulges show " cores" , whether bulges 
show nuclear sources. 

Bulges of disk galaxies have been targeted less often 
than ellipticals by the HST Peletier et al. (1999, here- 
after Paper I) analyzed a sample of 19 field SO-Sbc galax- 
ies using WFPC2 F450W, F814W and NICMOS F160W 
images, with the goal of obtaining bulge stellar popula- 
tion diagnostics. The combination of blue and NIR col- 
ors allowed them to put tight limits on the ages of bulge 
populations. Ages of SO to Sb bulges were found to be 
comparable to those of ellipticals in the Coma cluster, 
with a small age spread < 2 Gyr (Sbc bulges showed col- 
ors corresponding to younger ages). Nuclei were found 
to be dusty, with Ay — 0.6 — 1.0 mag. 

Carollo and collaborators surveyed mid- to late-type 
bulges using WFPC2 and NICMOS (e.g. Carollo 1999; 
Carollo & Stiavelli 1998; Carollo et al. 1997, 1998, 2001, 
2002; Seigar et al. 2002). These authors focus on bulge 
structure. They provide fits using the 14, exponential 
and Nuker models, and propose a structural classification 
of bulges into '14 classical' and 'exponential'. Carollo et 
al. (2002) find nuclear resolved components (NC) in the 
centers of 60% of the exponential bulges. In their view, 
'14' and 'exponential' bulges respectively show 'high' and 



'low' nuclear profile slopes, a structural difference which 
would trace different formation histories. 

Whether bulges come in two families with distinct 
structural properties has implications for formation 
mechanisms of bulges. Several models have been pro- 
posed (see Wyse, Gilmore, & Franx 1997; Bouwens, 
Cayon, & Silk 1999; Kormendy & Kennicutt 2004): early 
collapse (Renzini 1999; Zoccali et al. 2003); mergers prior 
to disk formation (Kauffmann, Chariot & White 1996); 
satellite accretion (Pfenniger 1993; Aguerri et al. 2001); 
and disk instabilities (Pfenniger & Norman 1990; Zhang 
1999). Bulges with 14 structure fit in the early collapse or 
merger scenarios, while exponential bulges are destroyed 
by mergers (Aguerri et al. 2001) and may instead be 
expected from disk instabilities (Combes et al. 1990). 
Edge-on, peanut-shaped bulges are known to have bar 
dynamics, and are therefore also expected to form from 
disk instabilities (Kuijken & Merrifield 1995; Bureau & 
Freeman 1999). The existence of two classes of bulges 
is commonly understood as evidence that massive bulges 
come from mergers while less massive bulges grow as a 
result of disk instabilities (see, e.g., Athanassoula 2005). 

In this paper we analyze the structural properties of 
bulges of early- to intermediate-type galaxies at HST 
resolution using the SO-Sbc sample presented in Paper 
I. We address profile shapes, nuclear sources, nuclear 
slopes, and central massive black hole mass estimates. 
Given the ability of the Sersic model to describe the pro- 
files of spheroids at ground-based resolution, we use the 
Sersic model as our starting point and enquire whether 
the increased spatial resolution of the HST contributes 
to support or to modify the ground-based picture. We 
perform a bulge-disk decomposition of the surface bright- 
ness profiles using combined TKT+ground-based profiles 
that sample the galaxy light distribution from the nu- 
cleus to the disk-dominated region. Ignoring this step 
would bring up two problems: the un- modeled disk con- 
tribution to the inner profile would bias the bulge nu- 
clear parameters; and, we would not be able to derive 
basic bulge parameters such as the total luminosity and 
the effective radius as the HST images do not cover the 
entire bulge at the distances of our target galaxies. 

We avoid using the 14 or exponential models, rather we 
focus on Sersic fits to the bulge profiles to test if the pro- 
file shape dichotomy appears when it is not forced. Our 
first results on bulge profile shapes using HST data were 
presented in Balcells et al. (2003, hereafter Paper II). In 
that paper we show that 14 bulge profiles are exceed- 
ingly rare. In this and a companion paper (Balcells, 
Graham, & Peletier 2007, hereafter Paper IV) we per- 
form a comprehensive analysis of those profiles. We will 
show that inner surface brightness profiles show excesses, 
over the best-fit bulge Sersic model, which can be suc- 
cessfully modeled by adding central unresolved sources 
and/or inner exponential components to the fitting func- 
tion (§[3]). Section d] provides details on the estimation of 
parameter errors through fits to simulated profiles. The 
subsequent sections analyze the properties of the nuclear 
excess light. Sect. 15.11 shows that the galaxies with ex- 
tended nuclear components closely match those with nu- 
clear disky isophotes, which suggests that the excess light 
in the surface brightness profiles comes from flattened 
components such as disks, rings or inner bars. Sect. 15.21 
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derives luminosities and masses for the unresolved nu- 
clear sources and addresses the Compact Massive Ob- 
ject (CMOs) paradigm, i.e., that nuclear star clusters 
are the low-mass extension to central supermassive black 
holes. In § 15.31 we relate the point sources to black hole 
mass estimates from the bulge velocity dispersions. Fi- 
nally, in § 15.41 we present and discuss the nuclear surface 
brightness profile slopes and compare them to those of 
ellipticals, bulges and dwarf ellipticals. In Paper IV, we 
discuss global bulge and disk scaling relations as inferred 
from the profile decompositions. A Hubble constant of 
Hq = 75 kms -1 Mpc -1 is used throughout. 

2. GALAXY SAMPLE AND DATA 

We have analyzed 19 galaxies from the Balcells & 
Peletier (1994, hereafter BP94) diameter-limited sam- 
ple of inclined, early-to-intermediate type disk galaxies. 
The BP94 sample was selected from the Uppsala Gen- 
eral Catalog of Galaxies (UGC; Nilson 1973) to include 
all disk galaxies of types SO to Sbc, listed as unbarred 
in the UGC, with blue diameters greater than 2 arcmin, 
inclinations above 50° (i.e. mid- to high-inclination), and 
apparent blue magnitudes brighter than 14.0 mag, within 
given limits of equatorial and Galactic coordinates. Upon 
inspection, some cases were excluded due to being obvi- 
ously barred, interacting, or very dusty, leaving 30 galax- 
ies that were analyzed in BP94, APB95, and Peletier & 
Balcells (1996, 1997). The present subsample comprises 
19 galaxies of types SO to Sbc that were imaged with HST 
(NICMOS-F160W [camera 2] and WFPC2-F450W and 
F814W; Paper I). The subsample was selected to pro- 
vide representative examples of each Hubble type, and 
to exclude cases where dust obscured the nuclei. Due to 
their high inclinations, some of the galaxies may harbor 
bars which go undetected in the images; bars may be 
suspected from the peanut-shaped isophotes of some of 
the bulges (e.g. Kuijken & Merrifield 1995). None of the 
galaxies have a Seyfert or a starburst nucleus. 

The sample has been extensively studied by us in pre- 
vious papers. Nuclear colors at HST resolution and 
bulge ages have been presented in Paper I, where postage- 
stamp images of the HST data, for the 19 galaxies studied 
here may be found. Peletier & Balcells (1997) published 
X-band surface brightness profiles and isophotal parame- 
ters from ellipse fits to wider field-of-view UKIRT images. 
Central stellar velocity dispersions and a Fundamental 
Plane analysis are given in Falcon-Barroso, Peletier & 
Balcells (2002). Minor-axis kinematic profiles are pre- 
sented in Falcon-Barroso et al. (2003). 

The inclined viewing angle for this sample presents 
advantages and disadvantages for a structural study of 
the nuclear properties of disk galaxies. The main draw- 
back is extinction, which can completely hamper detec- 
tion of nuclear structures at visible wavelengths. Fortu- 
nately, extinction is smaller at NIR wavelengths. Paper I 
concludes that extinction for this sample is on average 
Ah = 0.1 — 0.2 mag in the nucleus, and much lower fur- 
ther out. Problems related to extinction should therefore 
be minor for the study presented here. The main advan- 
tage of working with an inclined sample is that isophotes 
provide information on the flatness of each galaxian com- 
ponent, and hence, they guide in the identification of 
these components. 

Basic properties of the sample are given in Table [TJ 



where we list distances, spatial scales, K- and i?-band 
absolute magnitudes, central velocity dispersions, and 
disk ellipticities. Distances are derived from recession 
velocities relative to the Galactic-standard of rest, from 
de Vaucouleurs et al. (1991, hereafter RC3). We present 
new K-hanA apparent and absolute magnitudes for the 
program galaxies. The photometry we had published in 
APB95 is about 0.5 mag brighter than that from 2MASS; 
new ellipse fits to the APB95 images yield total if-band 
apparent magnitudes that closely match 2MASS total 
magnitudes (mean difference of 0.00 mag, rms of 0.21 
mag). The coincidence of our photometry and that of 
2MASS makes us believe that the present photomet- 
ric zero points are more accurate than those of APB95. 
Therefore we adopt the new apparent magnitudes in this 
paper, except for NGC 5746 and 5965, which overfill our 
frames, and for NGC 5879, whose frame suffers from 
bad sky planarity. For those three galaxies we adopt 
the total absolute magnitudes from NED. Adopted ap- 
parent magnitudes are listed in Table [1] Absolute K- 
band magnitudes are corrected for Galactic extinction, 
cosmological dimming, and K-correction. Errors in the 
-ftf-band absolute magnitude include the photometric er- 
ror and a distance modulus error which assumes a flat 50 
kms -1 recession velocity error. We use the disk elliptic- 
ities from APB95, which have been derived on i^-band 
images; these ellipticities show minimal differences, of at 
most 0.05, with respect to the i?-band values given in 
BP94. 

Details of the HST observations and the reduction 
of the HST data are given in Paper I. Here we de- 
rive elliptically-averaged surface brightness profiles and 
isophotal shapes from the HST/NICMOS F160W images 
(19"xl9", 0.075 arcsec/pixel), from 0.03 arcsec (1/2 a 
pixel size) to typically 10 arcsec, using the galphot pack- 
age (J0rgensen, Franx, & Kjaergaard 1992). We keep the 
centers fixed at the galaxy luminosity peak, and let the 
ellipticity and position angle of the isophotes vary. 

To extend the surface brightness profiles to large radii, 
we use the elliptically-averaged if-band surface bright- 
ness profiles derived from UKIRT /IRC AM3 images (mo- 
saics of 75"x75" frames, 0.291"/pbcel) published by 
Peletier & Balcells (1997), which we transform to ii-band 
by approximating the H — K profiles with the transfor- 
mation 

H-K = 0.111(7 -K)- 0.0339 (2) 

derived from the GISSEL96 models of Bruzual & Char- 
lot (see Leitherer et al. 1996) using I — K profiles from 
Peletier & Balcells (1997). 

The HST and ground-based (GB) profiles have match- 
ing slopes in the range 3" < r < 8" and show zero- 
point offsets that are always below 0.1 mag. We correct 
these by applying an offset to the GB profiles, which 
overall have less photometric accuracy. The process de- 
scribed here is the same that was employed in Paper I 
to derive color profiles for bulges at HST resolution. 
The match of GB to HST profiles is extremely good, 
as can be seen in the residual profiles presented below 
(see Appendix, Fig. lAl[) . The resulting geometric- mean- 
axis profiles were used in Paper II for the analysis of the 
Sersic shape index n. 

3. PROFILE DECOMPOSITION 
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TABLE 1 
Galaxy sample. 



NGC 


T 


D 


scale 


K 


± 


M K 


± 


Mr 


Co 


± 


EDisk 






Mpc kpcarcsec 


mag 


ma 


g 


mag 


km s 


-i 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) (11) 


(12) 


5326 


1 


34.3 


0.166 


8.99 


0.12 


-23.72 


0.13 


-21 


.28 


164 


6 


0.55 


5389 





26.2 


0.127 


8.57 


0.06 


-23.55 


0.08 


-20 


.11 


114 


6 


0.75 


5422 


-2 


25.3 


0.123 


8.88 


0.05 


-23.16 


0.08 


-21 


.08 


160 


6 


0.80 


5443 


3 


27.2 


0.132 


9.06 


0.23 


-23.13 


0.24 


-20 


.65 


76 


8 


0.72 


5475 


1 


24.5 


0.119 


9.30 


0.03 


-22.67 


0.07 


-20 


.10 


91 


6 


0.71 


5577 


4 


19.6 


0.095 


9.53 


0.08 


-21.96 


0.11 


-19 


.63 






0.72 


5587 





31.0 


0.150 


9.61 


0.09 


-22.88 


0.10 


-20 


.19 


93 


8 


0.70 


5689 





30.3 


0.147 


8.50 


0.28 


-23.94 


0.28 


-21 


.41 


113 


6 


0.75 


5707 


2 


31.1 


0.151 


9.31 


0.08 


-23.18 


0.09 


-20 


.55 


111 


6 


0.75 


5719 


2 


23.1 


0.112 


8.40 


0.09 


-23.44 


0.11 


-20 


.72 


108 


6 


0.68 


5746 


3 


22.9 


0.111 


6.88 


0.04 


-24.95 


0.07 


-21 


.80 


139 


8 


0.83 


5838 


-3 


18.3 


0.089 


7.72 


0.07 


-23.63 


0.11 


-21 


.01 


255 


6 


0.63 


5854 


-1 


23.3 


0.113 


8.63 


0.17 


-23.24 


0.18 


-20 


.55 


97 


6 


0.70 


5879 


4 


13.8 


0.067 


8.79 


0.22 


-21.92 


0.24 


-19 


,12 


58 


8 


0.70 


5965 


3 


47.6 


0.231 


8.61 


0.20 


-24.82 


0.20 


-22 


.04 


162 


8 


0.83 


6010 





26.0 


0.126 


8.82 


0.34 


-23.31 


0.34 


-20 


.69 


111 


6 


0.77 


6504 


2 


61.9 


0.300 


9.21 


0.07 


-24.81 


0.07 


-23 


.69 


185 


6 


0.80 


7457 


-3 


13.7 


0.066 


8.70 


0.16 


-22.00 


0.19 


-20 


.03 


56 


6 


0.48 


7537 


1 


35.3 


0.171 


9.68 


0.20 


-23.11 


0.20 


-20 


,19 


42 


9 


0.66 



Note. — Column description: (1) Galaxy NGC number. (2) Morphological type index from the RC3 (de Vaucouleurs et al. 1992). (3) 
Galaxy distance, computed from the Galactic-standard-of-rest recession velocities listed in the RC3, assuming H = 75 km s" 1 Mpc" 1 . (4) 
Spatial scale in kpcarcsec -1 at the galaxy distance. (5,6) Galaxy A"-band apparent magnitude and error, from our photometry (§[2j. (7,8) 
Galaxy K-band absolute magnitude, using the distances in column (2), with Galactic extinction, cosmological correction and K-correction. 
(9) Galaxy i?-band absolute magnitude, from BP94, for the assumed cosmology. (10,11) Aperture-corrected central velocity dispersion and 
error, from Falcon-Barroso et al. (2002). (9) Disk ellipticity derived from K-band images, from APB95. 



TABLE 2 

Best-fit parameters for the disk, bulge, and nuclear components. 



NGC 


A»o 


h 




Rc 




n 


B/D 


Hps 


H E2 


MO, 2 




h 2 


Fit type 


(1) 


(2) 


(3) 


(4) 


(5) 




(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


5326 


17.10 


11.47 


15.93 


3.88 


2 


.60 


0.99 




13.64 


11.55 





.15 


esc 


5389 


16.18 


9.70 


14.90 


2.29 


1 


.35 


0.39 




15.19 


11.04 





.06 


esc 


5422 


16.11 


8.61 


15.34 


2.56 


2 


.04 


0.47 


17.79 










pse 


5443 


16.67 


11.56 


15.28 


1.41 


1 


.69 


0.13 




16.34 


10.57 





.03 


ese 


5475 


16.15 


7.27 


16.32 


2.52 


2 


.12 


0.27 


17.51 


13.50 


12.24 





.22 


pese 


5577 


17.50 


13.30 


18.45 


2.75 


1 


.09 


0.04 












se 


5587 


16.28 


7.63 


15.69 


1.43 





.71 


0.10 




14.49 


13.59 





.26 


ese 


5689 


15.33 


7.38 


14.54 


2.03 





.86 


0.28 




13.39 


12.48 





.26 


esc 


5707 


16.25 


7.41 


15.53 


2.72 





.78 


0.45 


16.87 


12.95 


12.10 





.27 


pese 


5719 


17.30 


14.83 


15.72 


5.48 


1 


.96 


1.50 




15.34 


10.06 





.01 


ese 


5746 


15.24 


14.63 


15.15 


2.98 


1 


.55 


0.10 


17.65 










pse 


5838 


16.54 


17.24 


15.14 


4.82 


1 


.11 


0.63 


15.76 


11.42 


11.21 





.36 


pese 


5854 


16.15 


9.30 


16.44 


3.94 


1 


.88 


0.35 


15.45 


12.96 


11.02 





.16 


pese 


5879 


15.08 


5.48 


16.54 


2.18 


2 


.23 


0.11 


17.97 










pse 


5965 


16.21 


10.99 


16.08 


2.68 


2 


.75 


0.20 












se 


6010 


15.81 


8.08 


15.02 


2.04 


2 


.12 


0.36 


16.89 










pse 


6504 


16.22 


7.55 


16.32 


3.54 


2 


.65 


0.59 


17.85 










pse 


7457 


16.61 


14.50 


16.45 


3.41 


1 


.99 


0.17 


15.53 


14.55 


12.53 





.16 


pese 


7537 


15.64 


4.44 


17.14 


1.12 


1 


.76 


0.04 


17.99 










pse 



Note. — Table lists output parameters from the surface brightnesss profile fitting code, prior to applying the corrections described in 
§|4] All surface brightnesses are given in //-band magarcsec - 2 , and are not corrected to face-on view. All scale-lengths, in arcsec, refer to 
the geometric mean axis \f(a X b) of each measured ellipse. (1) Galaxy NGC number. (2) Extrapolated disk central surface brightness. (3) 
Disk scale length. (4) Bulge effective surface brightness. (5) Bulge effective radius. (6) Bulge Sersic index. (7) Luminosity ratio between 
bulge and main disk, from best-fit parameters (//-band). (8) //-band magnitude of central unresolved source, from best-fit parameters 
(//-band). (9) //-band magnitude of nuclear exponential component, from best-fit parameters (//-band). (10) Extrapolated central surface 
brightness of nuclear exponential. (11) Scale length of nuclear exponential. (12) Fit type code, se: Sersic bulge and exponential outer 
disk, pse: se plus a nuclear point source, ese: se plus an inner exponential component, pese: se plus a nuclear point source and a nuclear 
exponential component. 
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TABLE 3 

Physical parameters for the disk and bulge components. 



NGC Mx.Buigo ± M K , Disk ± B/D no ± log(h) ± ± log(i? c ) ± log(n) ± 

(1) (2) (3) (4) (5) (6) (7) (8) (9) (10) (11) (12) (13) (14) (15) (16) 



5326 


-22. 


.96 


0. 


17 


-22. 


.98 


0. 


1 1 


0. 


98 


17. 


94 


0. 


.04 





.455 





.004 


15. 


.69 


0. 


35 


-0. 


.23 


0. 


07 


0. 


11 


0. 


.07 


5389 


-22. 


.17 


0. 


14 


-23. 


.19 


0. 


.09 


0. 


39 


17. 


66 


0. 


.04 





.392 





.004 


11. 


66 


0. 


35 


-0. 


.58 


0. 


07 


0. 


13 


0. 


.07 


5422 


-21. 


.86 


0. 


10 


-22. 


.77 


0. 


.13 


0. 


.43 


17. 


84 


0. 


12 





.374 





.014 


15. 


.26 


0. 


22 


-0 


.53 


0. 


05 


0. 


31 


0. 


.01 


5443 


-20 


.76 


0. 


26 


-23. 


.00 


0. 


24 


0. 


13 


18. 


03 


0. 


.04 





.460 





.004 


15. 


,06 


0. 


35 


-0 


.77 


0. 


07 


0. 


23 





.07 


5475 


-20 


.98 


0. 


14 


-22. 


.41 


0. 


08 


0. 


27 


17. 


47 


0. 


.04 





.206 





.004 


16. 


10 


0. 


35 


-0. 


.57 


0. 


07 


0. 


33 





.07 


5577 


-18. 


.20 


0. 


13 


-21. 


.93 


0. 


15 


0. 


03 


18. 


85 


0. 


12 





.378 





.014 


18. 


.36 


0. 


22 


-0 


.61 


0. 


05 


0. 


01 


0. 


.04 


5587 


-20 


.26 


0. 


16 


-22. 


.78 


0. 


11 


0. 


10 


17. 


56 


0. 


.04 





.321 





.004 


15. 


.46 


0. 


35 


-0 


.71 


0. 


07 


-0. 


.15 


0. 


.07 


5689 


-22. 


.27 


0. 


31 


-23. 


.68 


0. 


29 


0. 


27 


16. 


80 


0. 


.04 





.337 





.004 


14. 


,31 


0. 


35 


-0 


.57 


0. 


07 


-0. 


.07 





.07 


5707 


-21 


.90 


0. 


15 


-22. 


.78 


0. 


10 


0. 


44 


17. 


73 


0. 


.04 





.350 





.004 


15. 


,30 


0. 


35 


-0 


,13 


0. 


07 


-0. 


.11 





.07 


5719 


-22. 


.88 


0. 


16 


-22. 


.15 


0. 


.12 


1. 


.49 


18. 


51 


0. 


04 





.467 





.004 


15. 


,48 


0. 


.35 


-0 


.26 


0. 


07 


0. 


29 


0. 


.07 


5746 


-22. 


.31 


0. 


10 


-24. 


.85 


0. 


.13 


0. 


10 


17. 


13 


0. 


12 





.596 





.014 


15. 


06 


0. 


22 


-0. 


.51 


0. 


05 


0. 


.19 


0. 


.01 


5838 


-22 


.59 


0. 


16 


-23. 


.10 


0. 


12 


0. 


62 


17. 


59 


0. 


.04 





.401 





.004 


14. 


,91 


0. 


35 


-0. 


,11 


0. 


07 


0. 


16 





.07 


5854 


-21. 


.75 


0. 


22 


-22. 


.92 


0. 


19 


0. 


34 


17. 


42 


0. 


.04 





.282 





.004 


16. 


20 


0. 


35 


-0. 


.40 


0. 


07 


0. 


27 





.07 


5879 


-19. 


.35 


0. 


25 


-21. 


.81 


0. 


26 


0. 


10 


16. 


38 


0. 


.12 


-0 


.174 





.014 


16. 


,47 


0. 


22 


-0. 


.87 


0. 


05 


0. 


35 





.04 


5965 


-22. 


.81 


0. 


21 


-24. 


.63 


0. 


23 


0. 


.19 


18. 


10 


0. 


12 





.789 





.014 


15. 


,98 


0. 


22 


-0. 


.24 


0. 


05 


0. 


14 


0. 


.04 


6010 


-21 


.79 


0. 


35 


-23. 


.00 


0. 


36 


0. 


33 


17. 


35 


0. 


12 





.328 





.014 


14. 


.90 


0. 


22 


-0. 


.62 


0. 


05 


0. 


33 


0. 


.04 


6504 


-23. 


.68 


0. 


10 


-24. 


.34 


0. 


13 


0. 


55 


17. 


91 


0. 


12 





.704 





.014 


16. 


20 


0. 


22 


-0 


.00 


0. 


05 


0. 


42 





.04 


7457 


-19. 


.87 


0. 


23 


-21. 


.84 


0. 


20 


0. 


16 


17. 


29 


0. 


.04 





.125 





.004 


16. 


22 


0. 


35 


-0 


.69 


0. 


07 


0. 


30 





.07 


7537 


-19. 


.47 


0. 


22 


-23. 


.07 


0. 


23 


0. 


01 


16. 


76 


0. 


12 





.114 





.014 


17. 


.03 


0. 


22 


-0. 


.75 


0. 


05 


0. 


21 


0. 


.04 



Note. — Table lists parameters for disk and bulge corrected from measurement offsets as described in §[4] Absolute magnitudes are 
given in the K-band. Surface brightnesses are given in H-band magarcsec - 2 . Both magnitudes and surface brightnesses are corrected 
for Galactic extinction, cosmological dimming, and K-correction. Scale-lengths are in kpc; for bulges, they refer to the geometric mean 
axis (ab) 1 ' 2 of the measured ellipse, while, for disks, scale-lengths are scaled to the major axis, assuming an inclination given by the disk 
cllipticity. (1) Galaxy NGC number. (2-5) Bulge and disk A*-band absolute magnitudes and errors, from the galaxy ii"-band absolute 
magnitude and the bulge-disk ratio from col. (6). (6) Luminosity ratio between bulge and main disk, from best-fit parameters (H-band). 
(7,8) Disk face-on extrapolated central surface brightness, and error. (9,10) Disk scale length, and error. (11,12) Bulge effective surface 
brightness, and error. (13,14) Bulge effective radius, and error. (15,16) Bulge Sersic index, and error. 



TABLE 4 

Physical parameters for the nuclear components. 



NGC M K ,ps 


± 


M K ,E2 


± 


W),2 


± 


log(h 2 ) 


± 


Fit type 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


5326 






-19.21 


0.47 


11.58 


0.94 


-1.60 


0.18 


ese 


5389 






-17.07 


0.47 


11.71 


0.94 


-2.00 


0.18 


esc 


5422 


-14.79 


0.35 














pse 


5443 






-15.99 


0.47 


11.12 


0.94 


-2.34 


0.18 


ese 


5475 


-15.00 


0.35 


-18.70 


0.26 


13.11 


0.45 


-1.40 


0.04 


pese 


5577 


















se 


5587 






-18.23 


0.26 


14.42 


0.45 


-1.23 


0.04 


ese 


5689 






-19.28 


0.26 


13.50 


0.45 


-1.20 


0.04 


ese 


5707 


-16.16 


0.35 


-19.77 


0.26 


13.13 


0.45 


-1.18 


0.04 


pese 


5719 






-16.64 


0.47 


10.45 


0.94 


-2.34 


0.18 


esc 


5746 


-14.72 


0.35 














pse 


5838 


-16.13 


0.35 


-20.16 


0.26 


11.81 


0.45 


-1.36 


0.04 


pese 


5854 


-16.96 


0.35 


-19.14 


0.26 


11.84 


0.45 


-1.56 


0.04 


pese 


5879 


-13.28 


0.35 














pse 


5965 


















se 


6010 


-15.78 


0.35 














pse 


6504 


-16.71 


0.35 














pse 


7457 


-15.71 


0.35 


-16.29 


0.47 


12.40 


0.94 


-2.02 


0.18 


pese 


7537 


-15.34 


0.35 














pse 



Note. — Table lists parameters for nuclear components corrected from measurement offsets as described in §|4] Nuclear source absolute 
magnitudes are given in the i\-band, using H — K = 0.23. Surface brightness fi0,2 is given in //-band magarcsec" 2 . Both magnitudes 
and surface brightnesses are corrected for Galactic extinction, cosmological dimming, and K-correction. Scale-lengths are in kpc, and have 
been scaled to the major axis, assuming an inclination given by the disk ellipticity. (1) Galaxy NGC number. (2,3) i\-band absolute 
magnitude of the central unresolved source, and error. (4,5) i\-band absolute magnitude of the nuclear exponential component, and error. 
(6,7) Face-on extrapolated central surface brightness, and error, of the nuclear exponential component. (8,9) Scale length, and error, of the 
nuclear exponential component. (10) Fit type code, se: Sersic bulge and exponential outer disk, pse: se plus a nuclear point source, ese: 
se plus an inner exponential component, pese: se plus a nuclear point source and a nuclear exponential component. 
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The decomposition of the surface brightness pro- 
files was described in Paper II. Briefly, the combined 
HST+GB profiles were fitted with a PSF-convolved 
Sersic plus exponential law using the code described in 
G01, modified to use a Moffat PSF {(3 = 6.9). The Moffat 
FWHM=0.131" used here is slightly narrower than that 
used in Paper II. This has two consequences: first the 
light in the NICMOS PSF wings needs to be accounted 
for with a suitable correction (§[3]); second, more nuclear 
components appear as resolved in the present fits than 
they did in Paper II. Fitting proceeded by chi-square 
minimization, allowing all 5 free parameters (disk [iQ,h; 
bulge fi e , R c , n) to vary. Convergence was generally 
straightforward, except for NGC 5326 and NGC 5854, 
where we fixed the disk parameters by eye. In a few 
other cases, a few outer los S/N points were excluded to 
prevent downward deviations from distorting the fit in 
the main part of the disk (see Figure[l|. 

We found that pure Sersic -(-exponential fits to the en- 
tire radial range of the profiles provide an inaccurate de- 
scription of the data: the residuals from the fits show 
a strong wave pattern with an obvious central positive 
residual. An example of this feature is shown in Figure 1 
of Paper II. As a consequence, the bulge and disk fitted 
parameters are highly sensitive to any inner radius cutoff 
imposed on the fitting range. Such fits yield uncertain 
values for the bulge Sersic index n, the total bulge lu- 
minosities, and the bulge-to-disk ratios. Fits excluding 
the inner ~0. 5-1.0 arcsec in radius (a common approach 
for deriving bulge parameters in the presence of nuclear 
components, e.g., Carollo et al. 1998; Stiavelli et al. 2001) 
show strong central residuals (see Fig. lAll second row). 
These indicate that our galaxy surface brightness pro- 
files cannot be modeled with just the sum of a Sersic 
bulge and an exponential disk. For 14 out of 19 galax- 
ies, we find an excess of central light above the Sersic 
bulge; 2 galaxies show central depressions, while 3 follow 
the Sersic profile reasonably well over the entire radial 
range. These numbers do not substantially vary when 
modifying the inner radial cutoff. 

The failure of pure Sersic -(-exponential fits suggests 
that many of our program galaxies harbor nuclear 
components. The detection of central excesses has 
previously been reported for bulges, dwarf ellipticals, 
and intermediate-luminosity ellipticals (Aaronson 1978; 
Binggeli, Sandage & Tarenghi 1984; Caldwell & Bothun 
1987; Phillips et al. 1996; Carollo et al. 1997; R01; Ravin- 
dranath et al. 2001). Excluding the central ~1" has been 
a common strategy to cope with these central compo- 
nents, which also avoids the problems with central dust 
extinction, e.g. the Stiavelli et al. (2001) fits to HST- 
based dwarf elliptical surface brightness profiles, or the 
14 or exponential fits to bulge surface brightness profiles 
by Carollo, Stiavelli, & Mack (1998). In our case, outside 
1", the pure Sersic -(-exponential bulge-disk fits are quite 
satisfactory (see Figure lATj) . suggesting that the Sersic 
model provides a good approximation to the large-scale 
brightness profiles of bulges. 

We therefore assume that the Sersic model describes 
"the bulge", taken as the spheroidal component residing 
in the center of the disk galaxy, and that any central 
excess above the Sersic profile is due to additional pho- 
tometric components. To implement our approach, we 



run the fitting program with additional central compo- 
nents. We tested central point sources, central exponen- 
tial disks, a combination of both, and central Gaussians 
of free width. In all cases, the models were convolved 
with the Moffat PSF prior to fitting. For each galaxy, 
all fits were inspected, and we selected the fits which, 
with a minimum number of added components, resulted 
in residual profiles without structures such as the wave 
pattern described above. In all cases, we adopted the 
extra component when the root-mean-square (rms) of 
the residuals improved by over 10%, and we discarded 
it when the rms improved by less than 10%. From a 
formal point of view, this approach is justified given 
that, with typically 90 independent data points, adding 
a point source to the 5-parameter Sersic+exponential fit 
would trivially yield an rms decrease of < 2.8%; includ- 
ing an additional nuclear exponential to the above so- 
lution would trivially improve the rms of that solution 
by 4.0%; and, an rms trivial improvement of 10% would 
be expected only after adding 18 constraints to the fit. 
Moreover, the adequacy of the employed model is ascer- 
tained through examination of the shape of the residual 
profile. Fits with residual profiles that were featureless 
were adopted without testing more complex models, as 
the rms cannot be improved by modifying the functional 
form of the fitted model in those cases. Residual pro- 
files for the adopted solutions have 0.02 < rms < 0.1 
magnitudes. 

Hence, our approach has two steps: first, a \ 2 
minimization of five different models, namely: pure 
Sersic+exponential; Sersic+exponential+Moffat; Ser- 
sic+exponential+inner exponential; Sersic+exponen- 
tial+Moffat+inner exponential; and Sersic+exponen- 
tial+Gaussian. Second, a selection of the best of these 
models, applying Occam's razor to choose the simplest 
model which describes the data. Our solutions describe 
the surface brightness profiles to a high degree of ac- 
curacy, although we generally cannot guarantee to have 
obtained unique solutions. Parameter uncertainties are 
discussed in § |U 

FigureQ] shows the profiles, the best-fit models, and 
the residual profiles from the fits for the 19 galaxies. Pa- 
rameters for the fits are given in Table [2] The largest 
structures in the residual profiles occur in the region of 
the disk, and are due to spiral arms and other disk fea- 
tures. In the bulge-dominated region, the small oscilla- 
tion around 0.2 arcsec arises from the cross-pattern in the 
NICMOS PSF wings and deviations from the assumed 
Moffat model (§gj. 

Internal extinction is obviously an issue at the galaxy 
centers. Fortunately, our use of NIR imaging data 
strongly alleviates the problem. Paper I shows that 
^4v,Contor = 0.6 — 1.0 mag on average for this sample, 
which corresponds to Ah between 0.1 and 0.2 mag. Ex- 
tinction is highly concentrated in the inner ^100 pc, 
hence is probably local rather than due to intervening 
disk dust. This extinction must vary from galaxy to 
galaxy, and can be estimated from the excess of central 
I — H color with respect to a typical bulge stellar popu- 
lation. The resulting values are so small that we do not 
apply extinction corrections to our photometry. The lack 
of such correction does not affect any of the conclusions 
of the paper. 
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Fig. 1. — H-band, combined HST plus ground-based surface brightness profiles for the 19 galaxies, and best-fit models. The abscissa is 
the geometric mean radius in arcsec, r = a — e), a = major axis, e = ellipticity. Filled circles: fitted data points. Open circles: outer 
points excluded from the fit. Solid lines: seeing-convolved Sersic bulge and exponential outer disk models. Dotted lines: seeing-convolved 
inner point sources. Dashed lines: inner exponential components. Given in each panel are the best-fit outer disk parameters (central 
surface brightness fio, disk scale length h) and bulge parameters (effective surface brightness fi e , effective radius R e and Sersic index n). 
Parameters for the inner components are given in Table [2] Below each profile is the profile residual (data minus model) ; A gives the 
root-mean-square residual from the fit, in magnitudes. 



4. PARAMETER UNCERTAINTIES 

Simulations with synthetic profiles were carried out to 
estimate errors in parameter recovery. The simulation 
work addresses two main questions. First is parameter 
coupling, especially between bulge and nuclear compo- 
nents, e.g., a profile consisting of a high-n Sersic and 
an outer exponential might be reproduced by the fitting 
program as a nuclear component, a lower- n Sersic and an 
outer exponential. A second important issue is whether 
our choice of a Moffat PSF is adequate. The Tiny-Tim 
PSF for our NICMOS images comprises a central peak 
surrounded by a secondary maximum at about 0.23". 
An analytical PSF such as the Moffat function is partic- 
ularly convenient for one-dimensional profile fitting, but 



may affect the parameters derived for nuclear sources, 
and the light in the secondary maximum may masquer- 
ade as an extended component. 

For the simulations, synthetic images were generated 
with IRAF's mkob j ects task, comprising a Sersic compo- 
nent, an outer exponential component, with and without 
nuclear components. Nuclear components were either a 
point source; a Gaussian; or an inner exponential. The 
sampled range for the parameters was bigger than that 
displayed by our target galaxies. The images were con- 
volved with Tiny-Tim PSFs derived from the HST NIC- 
MOS images of the target galaxies, and noise was added 
to yield surface brightess profile errors similar to those of 
the target galaxies. Surface brightness profiles were de- 
rived for each simulated galaxy image by fitting ellipses, 
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Fig. 1. 



continued. 



using the same fitting parameters that were employed 
for the derivation of the profiles for the program galax- 
ies. The radial extent of the profiles was somewhat lower 
than that of the program galaxies, but this should not 
affect the results given that, as we show below, disk pa- 
rameters were accurately reproduced with their current 
radial extent. The profiles were fitted with combinations 
of Sersic, exponential and nuclear components using the 
same code used for the program galaxies. As a test for 
any tendency of the fitting code to add non-existent nu- 
clear components, all models were fitted with and with- 
out nuclear components in the fitting function. 

We found that nuclear Gaussian components are par- 
ticularly difficult to reproduce. We suspect that the 
quadratic dependence of the Gaussian function on r 
made the fits unstable. Because of the failure with syn- 
thetic profiles, we do not present fits to the program 
galaxies employing Gaussian nuclear components. We 
note however that such fits yielded results that were gen- 



erally consistent with the results of PS or exponential 
nuclear components: galaxies well modeled with inner 
PSs yielded good fits with very narrow Gaussians, and 
galaxies that required an inner exponential component 
yielded good fits with broad Gaussians as well. 

Results from the simulations are summarized in Ta- 
ble El which lists mean offsets (measured minus input) 
and rms deviations for each group of models. For the 
statistics, we group together pure Sersic -(-exponential 
(se) models and PS+Sersic+exponential (pse) models, 
which show similar uncertainties, and list models with 
inner exponentials (ese) separately; differences between 
the two sets are small anyway. Disk parameters are re- 
covered with high accuracy. Uncertainties for bulge pa- 
rameters are only somewhat larger, and are asymmetric: 
we tend to recover fainter /x e (0.2-0.3 mag arcsec -2 ), and 
larger R c (A(log(i? c )) ~ 0.05). Such average uncertain- 
ties are probably overestimates, as the simulated models 
include inner components more luminous and more ex- 
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TABLE 5 

Parameter uncertainties 





SG ; pSG 


models 


ese models 


.r aram 


j\ 




A 




V 1 ) 






(4) 


(5) 


Outer exponential (disk): 










MO 


n n&o, 
u.uoy 


(119^ 
U. 1ZO 


-0.010 


0.036 


LOg^tl) 


u.uu / 




-0.003 


0.004 


Mag 


0.029 


0.085 


0.010 


0.046 


— 2 5 ]ae(D/T) 


0.066 


0.101 


0.012 


0.046 


Sersic (bulge) \ 












u.uuu 


224 


0.205 


0.348 


iog{n c ) 




n n^9 


0.044 


0.070 




-0.007 


0.035 


0.007 


0.072 


Mag 


-0.090 


106 


-0.019 


0.150 




-u. llO 


u.uy 


-0.025 


0.163 




-0.052 


0.071 


-0.017 


0.118 


Nuclear unresolved source: 










Magps 


0.310 


0.351 






Nuclear exponential, brigth range: 
















0.45 


0.45 


log(/i 2 ) 






0.09 


0.04 


Mag 2 






-0.003 


0.26 


Nuclear exponential, faint range: 










MO, 2 






0.81 


0.94 


log(/i 2 ) 






0.18 


0.18 


Mag 2 






-0.09 


0.47 



Note. — For each parameter, cols. (2) and (3) give the mean 
difference (measured minus input) and rms deviation for se, pse 
synthetic models. Columns (3) and (4) give the same quanti- 
ties for ese models. For the nuclear exponentials, bright range 
(faint range) denotes models in which the inner exponential is 
brighter(fainter) than the Sersic component at R = 0.1". See §|4] 



tended than seen in our program galaxies. 

Point sources are recovered with an offset of ~ 0.3 mag. 
Most likely, this offset arises from the fractional light in 
the secondary maximum of the NICMOS PSF (0.5 mag 
outside 2 pixels for F160W; Holfeltz & Calzetti 1999). 
The width of our Moffat PSF (FWHM = 0.131") was 
set to match the main peak of the NICMOS PSF, hence 
the light in the wings gets unaccounted for; a MOFFAT 
width of ~0.19" would greatly reduce such offset, but at 
the price of losing spatial resolution, therefore we adopt 
the narrower PSF and simply apply a 0.31 mag aperture 
correction to the point source magnitudes. 

For inner exponentials, models with ^0,2 >1 mag 
fainter than the Sersic's fto (Sersic) get lost in the noise, 
while, at the bright end, models with /io.2 — fJ-o (Sersic) < 
—2 yield inner profiles dominated by the inner exponen- 
tial, which are quite unlike those of real galaxies. The 
statistics shown in Table [S] correspond to the range be- 
tween those limits, which we split in bright and faint 
ranges as shown in the Notes to the table. The PSF wings 
affect inner exponential components by making their de- 
rived scale-lengths about 10-20% higher, while making 
the central surface brightness fainter. Total magnitudes 
show an uncertainty of a few tenths of magnitude - larger 
for the fainter components as expected. 

Three main lessons derive from the simulations, (i) 
Bulge and outer disk parameters are robust to the pres- 
ence of nuclear components and to the choice of PSF; in 
particular, our fitting code does not artificially introduce 



nuclear components to pure Sersic -(-exponential profiles 
with a high-n Sersic component, (ii) A moderate level 
of parameter coupling is present for the bulge: derived 
fi e are probably faint by ^0.2-0.3 magarcsec -2 , while 
log(i? ) are probably overestimated by ~5%. (Hi) Nu- 
clear unresolved sources are accurately recovered once an 
aperture correction of ~0.3 mag is applied. 

5. RESULTS 

Physical parameters for bulges, disks, and nuclear com- 
ponents are listed in Tables [3] and [4] Corrections for 
measurement offsets, derived from the simulations O, 
have been added where applicable (Table [5]). RMS error 
estimates derive as well from Table [5] 

The first result of our analysis is that nuclear sources 
are quite common in disk galaxy bulges. Excesses over 
the Sersic fit to the bulge are measured in 17 out of 19 
galaxies. Over half (11) of those are extended; six can 
be modeled with an exponential profile, while five harbor 
a central unresolved source in addition to the exponen- 
tial. The remaining 6 galaxies with central excesses can 
be modeled with an unresolved source. These results 
are not new. Nuclear sources, resolved or unresolved 
at the scale of the HST instruments have been reported 
by others, for bulges (Carollo et al. 1998; Paper I), for 
ellipticals and early- type bulges (R01; Ravindranath et 
al. 2001; F06), for dwarf ellipticals (Stiavelli et al. 2001; 
Graham & Guzman 2003; Cote et al. 2006), and for late- 
type spirals (Boker et al. 2002). The variety of detection 
methods and functions used by these teams to model the 
underlying light distributions suggests that nuclear com- 
ponents are not an artifact of the model fitting but are 
real components of many galaxy nuclei. 

Absolute magnitudes of the nuclei are plotted in Fig- 
ure [2] against if-band absolute magnitude of the bulge, 
bulge central velocity dispersion, bulge color and disk el- 
lipticity. Nuclei are on the mean 4.4 ±1.9 mag fainter 
than their host bulges, and 6.1 ± 1.7 mag fainter than 
their host galaxies, although some of the extended com- 
ponents are only 2 mag fainter than their host bulge (e.g., 
NGC 5587, see Table 2]). In general, nuclei are minor 
contributors to the galaxy light. 

If structural components are linked to formation events 
in the history of the host galaxies, we may enquire 
whether nuclei are recent additions to the galaxy, or are 
they old, perhaps the seed of the formation of the bulge or 
of the entire galaxy. If galaxy centers host a supermassive 
black hole, whose formation and growth went through 
a phase of positive feedback with the formation of the 
galaxy (Silk 2005), are the nuclei detected here connected 
in any way to such positive feedback (McLaughlin et al. 
2006)? In the following subsections we further analyze 
the extended and unresolved nuclear components. 

5.1. Nuclear extended components 

The nuclear extended components have absolute K- 
band magnitudes in the range —16 > Mj^Ext < —20. 
Their scale lengths are ~5 to 60 pc, or about 100 
times smaller than the outer disk scale lengths and 20 
times smaller than the effective radii of the bulge (for 
NGC 5719, an extreme case, i? e is 100 times higher 
than h^). But they have very high densities. Ex- 
trapolated central surface brightnesses are in the range 
11 < fin < 15, typically 5 mag brighter than the (j,q of 
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Fig. 2. — Absolute A"-band magnitude of the nuclear components against: (a) the A"-band absolute magnitude of the bulge. Filled 
circles: unresolved components in bulges. Open circles: resolved components in bulges. Triangles: nuclear components in dwarf ellipticals 
from Graham & Guzman (2003). The dotted line is an orthogonal regression to the bulge unresolved components, see eqn. [5] while the 
solid line is an orthogonal regression to the bulge and dE unresolved components together (eqn.[6]l. (b) aperture-corrected central velocity 
dispersion, from TableQ] (c) Central I-H color, from HST/NICMOS F160W and WFPC2 F814W images (Paper I). An extinction vector 
for normal Galactic extinction of Ay = 1 mag (Rieke & Lebofsky 1985) is plotted in the top-left, (d) Ellipticity of the outer disk, from 
two-dimensional bulge-disk decomposition in A"- band images (APB95). Upper limits are given for the galaxies without detected nuclear 
components; one of those (NGC 5577) does not have a central velocity dispersion measurement. 
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Fig. 3. — Disk scale-length in kpc vs. extrapolated central surface 
brightness, for inner extended components, and for main galaxy 
disks. Dotted line: orthogonal regression to the nuclear resolved 
components (Eqn. |3). Dashed line: orthogonal regression to both 
nuclear resolved components and outer disks (Eqn. l4t. 



the main galaxy disk. In many cases, nuclear extended 
components dominate the surface brightness profile in- 
ward of 1" (see Fig.HJ. 

On the basis of their luminosity distribution (Fig. [5^) 
the nuclear extended components might be a heteroge- 
neous family. Seven galaxies draw a bright sequence with 



absolute magnitudes — 18 > ^KKxt 

> -20, or 3% to 

15% of the bulge luminosity, while the four remaining 
cases show fainter luminosities (—16 > Mx Ext > -17, 
or 0.3% to 3% of the bulge luminosity) that overlap with 
the luminosities of the unresolved sources. These fainter 
objects are also the smallest {h^ < 10 pc) and densest, 
and may be instances of nuclear clusters that we man- 
aged to resolve in our images. The objects in the brighter 
sequence are also more extended (25 < ft. 2 < 60 pc), five 
of the seven cases themselves harbor nuclear unresolved 
sources, and their sizes and fractional luminosities re- 
semble those of inner disks in elliptical galaxies (Scorza 
& Bender 1995; Scorza & van den Bosch 1998). These 
properties give clues that such objects might have dif- 
ferent nature and formation mechanism than their more 
compact counterparts. The presence of more extended 
components than nuclear clusters indicates that struc- 
tural deviations from the Sersic profile in the nuclei of 
early-type disk galaxies are not restricted to nuclear star 
clusters: more extended structures, which we show be- 
low to be flattened systems, also cause the bulge profile 
to deviate from the Sersic functional form. 

Scale lengths and central surface brightnesses for inner 
and outer disks are plotted together in Figure [3l They 
correlate quite well, both for the nuclear disks alone, 



Io,2/Lk,q — 10 

{R s = 0.83; P nuU = 8 
disks together: 

I /L K , e = 10 



-3.12±0.16 



(/ia/10 pc) 



-1.26±0.31 



(3) 



9-10 ), and for inner and outer 



-5.30±0.12 



(fc/kpc) 



-1.07±0.05 



(4) 



(R s = 0.91; Pnuii = 8.4 • 10 ). In the above equations, 
Iq denotes central intensity in the if-band, in units of 
-fT-band solar luminosities per square arcsecond. Both 
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relations are consistent with each other, suggesting that 
Iq ~ h , or L ~ h. Caution is needed when interpret- 
ing these relations, given the important selection effects 
operating on the detection of inner components. We saw 
in |j4]that detection of inner disks is broadly constrained 
at —2 < /zo,2 — Mo,Sersic < 1- I n practice, we have found 
(Mo.2 — Mo,Scrsic) = —0.05 ± 0.81, i.e., only those inner 
disks that match the central surface brightness of the 
bulges's Sersic profile are detected, and any inner disks 
much fainter than that limit would go undetected by our 
fitting code. Hence, equations [3] and [4] describe the up- 
per envelope of inner disk surface densities. Similarly, 
the distribution of our large-scale disks in this diagram 
is known to define an upper envelope of points (Graham 
& de Blok 2001). 

We seek clues on the nature of nuclear extended 
components by relating their structural parameters to 
isophotal, dynamical and color information. Figures 
with such information are shown for each galaxy in Ap- 
pendixlAl 

All of the 11 galaxies with inner extended components 
show corresponding disky isophotes in the inner arcsec 
(see Fig. lAlj) . The association of extended components 
with disky isophotes is nearly one-to-one, as only three 
galaxies with disky isophotes in their bulges (NGC 5746, 
5965, and 7537) do not require an inner exponential com- 
ponent to fit the surface brightness profile. NGC 5746, 
with a peanut-shaped bulge, shows disky isophotes in the 
bulge-outer disk transition region; NGC 5965, another 
prototype peanut-shaped bulge, has very strong positive 
C4 disky coefficient in the entire bulge region, suggest- 
ing that the entire bulge is disk- like (see, e.g., Bureau 
et al. 2006); finally, NGC 7537 (Sbc), has a faint bulge 
and low velocity dispersion (a = 42 kms -1 ). Hence, the 
bulges in these three galaxies really look like disks, and 
fall under the definitions of pseudobulge (Kormendy & 
Kennicutt 2004). For the rest of the sample, inner disky 
isophotes are associated with nuclear extended compo- 
nents. 

Such correspondence suggests that inner extended 
components of the surface brightness profiles trace true 
structural components of the galaxy nuclei, with flat- 
tened shapes. On the basis of SAURON 3D spectroscopy, 
similar flat kinematically decoupled components aligned 
with the main galaxy disk have recently been reported 
for early-type spirals by Falcon-Barroso et al. (2006; see 
also McDermid et al. 2006). The simple conjecture that 
they are inner disks, bars, or nuclear rings, is supported 
by the kinematic data, which show velocities consistent 
with circular motions and often show velocity dispersion 
minima. Smaller, but also extended, 'flattened clusters' 
have been found in imaging programs for later type spi- 
rals by Seth et al. (2006). For our galaxies, the typical 
outer radius of the positive c4 structures is a few hun- 
dred pc, which is a typical size of inner bars in double-bar 
galaxies (Erwin & Sparke 2002, Erwin 2004). 

For our galaxies, nuclear features in the velocity dis- 
persion profiles (see Fig. lAlj) are not associated with the 
presence of nuclear extended components. In the inner 
arcsec where typically nuclear components have a strong 
contribution to the total surface brightness, most velocity 
dispersion profiles are quite flat, with specific instances 
of central peaks or drops. These profiles, from minor-axis 
spectra taken with the WHT/ISIS spectrograph, have a 



0.3 arcsec sampling and a seeing of around 1.0 arcsec 
(Falcon-Barroso et al. 2003). Clearly, two-dimensional 
spectroscopy at sub-arcsecond resolution is required for a 
kinematic characterization of the nuclear extended com- 
ponents. 

The colors of the extended nuclear components are 
quite red, and, as shown in Paper I, they correspond 
to a stellar population reddened with a mean of Ay — 
0.5 — 1.0. The color profiles, however, are generally very 
smooth (see Fig. lAip . This suggests that the stellar pop- 
ulations of the nuclear extended components are not very 
distinct from those of the surrounding bulges. In partic- 
ular, despite signs of nuclear star formation (Paper I), 
the bulk of their stellar populations are not recent addi- 
tions to the galaxy. Note that the combination of blue 
and NIR colors is particularly sensitive to population age. 
In later-type spirals, the nuclear 'flattened clusters' are 
bluer than the host spheroid (Seth et al. 2006), perhaps 
indicating an extended star formation history for the nu- 
clear structures which continues today for late-type spi- 
rals (see also Walcher et al. 2006), but has finished in 
most earlier type disk galaxies. 
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Fig. 4. — The spheroid mass is plotted against the mass of the 
compact massive object. Filled circles: this work, bulge mass vs 
nuclear unresolved source mass. Triangles: Coma dwarf ellipti- 
cals from Graham & Guzman (2003). The above are photometric 
masses, in Solar units, deriv ed fr om the _R"-band luminosities as- 
suming M/Lk = 0.8, see § 15.21 Open circles: ellipticals and 
bulges with central SMBHs, from Haring & Rix (2004). Solid line: 
regression for bulges and dEs (Eqn. 0. long-dash line: bisector 
linear regression for the sample of Haring & Rix (2004), see equa- 
tion [9] dotted lines: ±1<t range of acceptable fits. Dot-dashed 
line: the relation between CMOs and galaxy mass derived by F06 
under the assumption of a linear scaling (Eqn. [8}. 



5.2. Nuclear unresolved sources 

As previously noted, we detect unresolved sources (re- 
ferred to as point sources, or PSs) in 11 of our 19 galaxies; 
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five of these coexist with an inner exponential. Our de- 
tection frequency is similar to that reported in previous 
studies: Ravindranath et al. (2001) report a 50% detec- 
tion fraction; Boker et al. (2002) find nuclear PSs in 76% 
of their sample; Cote et al. (2006) report a fraction of 
66% to 88% of nucleated early-type galaxies from the 
ACS Virgo Cluster Survey. 

For our sample, the absence of AGN nuclei, and the 
nuclear colors suggest that light from the PSs is stellar 
in origin, i.e. that they are nuclear star clusters. A simi- 
lar conclusion is reached by Phillips et al. (1996) and by 
Carollo ct al. (1998), whereas Ravindranath et al. (2001) 
argue for a non-stellar origin for the PS light on the ba- 
sis of the high frequency of AGN in their sample. In 
what follows we assume the PSs to be nuclear clusters, 
although we recognize that, in general, AGN contribu- 
tion to nuclear emission may be important in specific 
samples. 

Absolute magnitudes are in the range — 13 > Mk, ps > 
— 17, (Table 2]), on the mean 6 mag fainter than their 
host bulges (if -band luminosity ratio 0.4%). Their lu- 
minosities correspond to 10-20 globular clusters, quite 
comparable to nuclei in nucleated Virgo ellipticals (Cote 
et al. 2006). These absolute magnitudes do not vary sys- 
tematically with galaxy inclination (Fig.[2ji), suggesting 
that extinction in the parent disks do not affect Mk,ps- 
Extinction at the nuclei themselves does have a small ef- 
fect on Mk, ps- We show for reference an Ay — 1 mag 
extinction vector in Figure^;. 

We showed in Paper II that Mk,ps correlate with the 
bulge absolute magnitude, a result later confirmed by 
Rossa et al. (2006) using optical imaging. The correla- 
tion given in Paper II, eqn. 1, was derived from fits to 
17 objects, using a slightly wider PSF (0.19" FWHM) 
than used here (0.131"). For the fits presented here, an 
orthogonal regression to the M^ps-M^buI distribution 
gives 

t it i n7.70i0.17/ t /in 10 r \0.63±0.37 

^K.PS/ ^K,Q = 1U (.-kK.Bul/J-U ^K,Q) , 

(5) 

(using Mk,q = 3.41, Allen 1973), consistent with that we 
gave in Paper II. However, with only 11 data points and a 
significant scatter, the relation is not statistically signif- 
icant for this sample (R s = 0.43; P nu ii = 0.18). Graham 
& Guzman (2003), who analyzed Coma dwarf elliptical 
(dE) i?ST/WFPC2 F606W galaxy surface brightness 
profiles using the same fitting code as the present pa- 
per, also find a scaling between PS and bulge luminosity, 
with a consistent slope of 0.87 ± 0.26. Their data points 
are given in Figured 5 . Inasmuch as the nuclei of our 
bulges and those of the Coma dEs are similar structures, 
we may use the combined sample to derive the scaling of 
nuclear unresolved sources with spheroid luminosity. We 
find 

t it nn 7.75±0.15/ r /in!0 T \0.76±0.13 

(6) 

(R a = 0.72; P uu u = 4.0 • 10~ 4 ). We transform eqn. ©into 
a mass relation, i.e., 

5 After applying a constant color term F606W — K = 2.7. This 
scaling corresponds to an old population with 0.4 times solar metal- 
licity, using the models of Vazdekis et al. (1996); a lower metallicity 
would make the total and nuclear luminosities of the dE fainter. 



Mp S /M & = 10 7 - 73±0 - 16 (A4bui/10 10 X q )°- 76±0 - 13 , (7) 

where M/Lk = 0.8 has been assumed for nuclei and 
bulges, from Bell & de Jong 2001, for the typical colors 
of bulge populations. The normalization is insensitive to 
the choice of M/Lk because its slope is close to unity; us- 
ing the more extreme M/Lk — 0.5 would yield a nearly 
identical normalization term of io 7 - 68±016 . 

Equation [7] may be compared to scaling regressions 
found in other studies of nuclear clusters. F06 find 
■Mnuc ~ A4 B ^f ±a25 , which is ~ 2a from our slope. And, 
for their so-called compact massive objects (CMOs), 
a class encompassing central supermassive black holes 
(SMBHs) and central star clusters, F06 find, when fixing 
the exponent to unity, 

Mcmo/M q = 10 7 - 26±a47 (X Gal /10 10 M 9 ) L0 . (8) 

while, for a sample of galaxies with kinematic SMBH 
mass determinations, Haring & Rix (2004, hereafter 
HR04) infer 

M./Me = 1O 7 - O8±o - 1o (Xbui/1O 1o X ) 1 - 12±0 - 06 . (9) 

Figure 2] shows the host spheroid mass against the nu- 
clear mass of unresolved bulge and dE components. For 
comparison we plot the HR04 distribution of A4, vs host 
mass. The regressions provided by HR04 (eqn. [5]), F06 
(eqn. [SJ, and us (eqn. [7|) are plotted as well. The dis- 
tribution of A^ps converges with that of the M, regres- 
sion at the high-mass end, but it progressively deviates 
as we move to lower spheroid masses. Both our distribu- 
tion of A4ps and HR04's distribution of M., show a large 
scatter at log(AlBui) < 10, hence the slopes of the regres- 
sions have large uncertainties in that mass domain. How- 
ever the nuclear cluster masses are systematically above 
the extrapolation of the HR04 or F06 relations. From 
equation[71 nuclei fractional mass (M.ps / Mbu\) increases 
from 0.19%, identical to central black holes fractional 
masses, to 0.71%, as A4b u i decreases from l0 ll b M.Q 
to 1O 9 5 A / 10. We are aware that different mass deter- 
minations are being compared in Figure [5] (photometric 
masses for our data, virial masses for F06, and masses 
derived from Jeans equations' modeling for HR04). Nev- 
ertheless, we argue that the observed offset between our 
unresolved components and the SMBH relations from ei- 
ther HR04 or F06, cannot be explained by uncertain- 
ties in our mass determinations. Bridging the offset by 
changing our assumed M/Lk would require M/Lk of 
nuclei to be lower than those of their parent galaxies 
by a factor of ~4, which is implausible for stellar pop- 
ulations without signs of vigorous star forming activity. 
Using dynamical masses for our bulges, instead of the 
photometric masses, does not help either: applying for 
consistency the same form of the Virial mass used by 
F06, 

A^BuLDyn =5i? <7 2 /G (10) 

the offset actually increases (A^Bui.Dyn are on the mean 
a factor 2 below the stellar-mass A^Bui) which shifts our 
points to the left, away from the F06 relation; such dis- 
crepancy indicates that eqn. 1101 underestimates the true 
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dynamical mass of the bulges - most likely due to ignor- 
ing the rotational kinetic energy. 

We thus cannot reproduce F06's result that nuclear 
star clusters masses fall onto the same linear relation as 
defined by more massive central black holes. Any CMO- 
bulge mass relation that encompasses both central black 
holes and nuclear star clusters must be non-linear for 
bulge masses of say .MbuI < 10 10 Mq. An orthogonal 
regression to the nuclear star cluster masses and the black 
hole masses, against the host bulge mass, gives 



M 



CMO 



/M c 



10 



7.51±0.06 



(AWlO^Afo 



\0.84±0.06 



(11) 

We conjecture that F06 find a linear A^cmo _ A^bu1 re- 
lation because, near the low mass end, their sample in- 
cludes many SO's. When deriving Virial masses using 
equation I10| because their R c derive from single Sersic 
fits to the galaxy profiles (Cote et al. 2006), spheroid 
masses are probably overestimated for all of their SO 
galaxies, which is roughly half of their sample. Lau- 
rikainen ct al. (2005) report that the mean bulge-to-total 
ratio for lenticular galaxies is 0.24, so one may expect 
a downward correction to (half of) the F06 spheroidal 
masses of ~0.6 dex, which would explain the different 
trends in F06 and us. 
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Fig. 5. — Central supermassive black hole mass, estimated from 
the velocity dispersion following Tremaine et al. (2002), plotted vs. 
(a) the photometric bulge mass, and (fe) the photometric mass of 
the nuclear unresolved sources. Photometric masses are der ived 
from the if-band luminosities assuming M/Lj^ = 0.8, see § 15.21 
Th e solid lines in panels (a) and (b) correspond to eons. [T3l and 
1141 respectively, while dotted lines delineate the approximate region 
of acceptable fits given the l-o~ error bars in the fit coefficients. In 
panel (a), the dot-dashed line traces the SMBH mass — shperoid 
mass relation of HR04, while the dashed line traces a slope unity 
relation for reference. In panel (b), the dashed line is the locus 
where PS masses and SMBH masses are equal. 



5.3. Nuclear black hole masses 



To further investigate the connection between star 
clusters and central SMBHs proposed by F06 and 
Wehner & Harris (2006), in this subsection we explore 
the scalings of SMBH mass A4 m , nuclear cluster mass, 
and bulge mass. We estimate A4, using the expression 

log M. = 4.02(±0.32) log(cr/200) + 8.13(±0.06) (12) 

from Tremaine et al. (2002 hereafter T02), where Ai, is 
in Solar masses, and a is the central velocity dispersion 
in kms -1 . We show bulge masses against SMBH masses 
in Figure EH- An orthogonal regression of our estimated 
M, against Mbm\ gives 

M./M Q = 1O 7 - 25±O - 11 (AW1O 1O ^ ) 1 - 72±O - 26 , (13) 

(i? s = 0.84; P null = 5.6 • 10~ 4 ). Such a strong rela- 
tion is expected because it is largely a manifestation 
of the strong Faber- Jackson relation followed by our 
bulges. Our bulges follow i^Bui oc (j 2 - 86±0 - 5 (see Pa- 
per IV, § 3.1.5.), and using M. oc cr 4 02±0 - 32 we ex- 
pect M. oc (A^bui) 4 02/2 - 86 = (-Mbui) 1 ' 41 ^- 31 , consis- 
tent with eqn. [T3] Equation [T5] is also consistent with 
that found by Laor (2001) for a sample of active and in- 
active galaxies, i.e., M., oc (Mbui) 1 ' 53 ^ ' 14 ■ Our scaling 
lies 2a from the relation found by HR04 (drawn in Fig- 
ure EH), and formally departs 2.8<7 from linearity (shown 
with a dashed line in Figure Eh.). We do not want to give 
much weight to this departure from linearity, given the 
small sample size, and the fact that we are extrapolating 
the A4 t — a relation faintward. 

If central black holes and central star clusters are in- 
timately related, the latter perhaps being failed black 
holes, then the A4,~ <r relation (eqn. [T2|) should be able to 
predict the nuclear cluster masses. We plot A^ps (mea- 
sured) against A4, (estimated) in Figure [SJd, together 
with the result of an orthogonal regression to the nuclear 
star-cluster-black hole mass distribution, which gives 

M PS /M Q = 1Q 7A6±0 - 14 (M./10 7 M G )°- 20±0 - 16 , (14) 

(R s — 0.37; Pnull = 0.24). The correlation is not statisti- 
cally significant, partly due to the small sample size and 
the scatter of the data points. We report this regression 
to clarify that the slope of the relation is several sigma 
away from the value of 1 expected if the A4,—o~ relation 
had predicted the masses of all the star clusters. This 
occurs also when SMBH masses are computed following 
Merritt & Ferrarese (2001). 

Wehner & Harris (2006) advocate a CMO transitional 
mass at 10 7 A4q. While our Ai, scatter above and below 
this mass limit, all but one of our nuclear star clusters lie 
above it. Some bulges are already known to contain both 
a SMBH and a nuclear star cluster: NGC 7457, which is 
in our sample, and NGC 3384 (Wehner & Harris 2006). 
Curiously, above the 1O 7 .M0 limit our nuclear star clus- 
ter masses agree with the predicted SMBH masses. Be- 
low the 10 7 7Wo limit, cluster masses lie ~1 dex above 
the extrapolation of the black-hole mass relation (shown 
as a dashed line in FigureEb). Hence, again, the M,-a 
relation appears not to predict the nuclear star cluster 
masses below this limit. It is unlikely that our surface 
brightness profile fitting code has overestimated the PS 
light by such large factors (§|3|, or that the M/L we 
adopted for the nuclear clusters is off by similar amounts. 
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Some of the deviation may arise from the large intrinsic 
scatter in the M, — a distribution in the low-mass range 
6 < \og(M.) < 8; see, e.g., Figure 7 from T02. 
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Fig. 6. — Inner surface brightness profiles plotted against log(i?) 
in the range 0.1 < R < 1 arcsec (R = geometric mean radius). Pro- 
files are sorted by decreasing bulge absolute magnitude, brightest 
above, and are offset by 1.5 mag/arcsec 2 from each other for clar- 
ity. Solid lines: direct least-squares power-law fits, with exponent 
7 given in Table[6] col. 2. 



5.4. Inner profile slopes 

The -ffST-resolved, logarithmic slopes of inner surface 
brightness profiles are a useful galaxy parameter as they 
can test the applicability of density profile models, which 
are often based on power-laws (e.g., Jaffe 1983; Hern- 
quist 1990; Dehnen 1993; Tremaine et al. 1994). Inner 
slopes constrain the shape of the potential and the types 
of orbits that may be present in the nuclei. On the 
observational side, interest in profile slopes arose from 
the discovery of a bimodal distribution of inner slopes 
in samples of intermediate- and high-luminosity ellipti- 
cal galaxies (Ferrarese et al. 1994; L95; Gebhardt et al. 
1996) which suggested different formation/evolution pro- 
cesses for the cores of giant and less luminous ellipticals, 
perhaps linked to the presence of binary black holes (e.g., 
F97). 




Fig. 7. — The spheroid B-band absolute magnitude plotted 
against two inner profile slope determinations 7. Top panel: slope 
7^ of the best-fit Nuker model at 0.3 arcsec. Bottom panel: slope 
7' of best-fit Sersic profile at 0.15 R e - Filled circles: bulges, this 
work. Squares: core ellipticals from L95. Triangles: power-law el- 
lipticals from L95. Five-pointed stars: dwarf ellipticals from Gra- 
ham & Guzman (2003). Six-pointed stars: dwarf ellipticals from 
Stiavelli et al. (2001). 



TABLE 6 
Inner negative logarithmic slopes. 



NGC 7(pwl) 


7(pwl) 




T S ' 


?8 




0.1-1" 


20-200 pc 


0.3" 


0.3" 


0.15R C 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


5326 


0.91 


0.85 


0.98 


0.90 


0.70 


5389 


0.69 


0.70 


0.72 


0.43 


0.39 


5422 


0.62 


0.65 


0.63 


0.72 


0.64 


5443 


0.83 


0.83 


0.84 


0.59 


0.72 


5475 


0.96 


1.03 


1.34 


0.75 


0.68 


5577 


0.17 


0.22 


0.16 


0.30 


0.22 


5587 


0.56 


0.52 


0.81 


0.11 


0.17 


5689 


0.46 


0.43 


0.49 


0.18 


0.17 


5707 


0.79 


0.77 


0.86 


0.14 


0.10 


5719 


0.66 


0.65 


0.67 


0.70 


0.42 


5746 


0.43 


0.46 


0.46 


0.53 


0.41 


5838 


0.67 


0.84 


0.78 


0.47 


0.26 


5854 


1.24 


1.26 


1.21 


0.67 


0.46 


5879 


0.63 


0.67 


0.65 


0.79 


0.76 


5965 


0.78 


0.61 


0.84 


0.94 


0.85 


6010 


0.74 


0.77 


0.76 


0.75 


0.75 


6504 


0.73 


0.54 


0.78 


0.92 


0.74 


7457 


0.96 


0.97 


0.92 


1.12 


0.82 


7537 


0.59 


0.63 


0.59 


0.62 


0.86 



Note. — Column description: (1) galaxy NGC number. (2) 7 
from power-law fit to 0.1"-l"(Fig.[6JI. (3) 7 from power-law fit to 
20pc-200pc. (4) 7' slope of Nuker model at R = 0.3" (single 
Nuker fit to 0.1" < R < 4"). (5) 7', slope of Sersic model at 
R = 0.3" (simultaneous bulge, disk and nuclear component fits). 
(6) 7', slope of Sersic model at R/R c = 0.15. 
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For our galaxies, all surface brightness profiles continue 
to rise inward to the resolution limit of the data (Fig.[6|). 
Inner profiles approximate power laws, although small 
but clear deviations are evident in Figure[6l making the 
slope determination dependent on the radial range used 
for the measurement. The presence of inner components 
further complicates matters and means that two types 
of profile slopes may be measured: the slope of the total 
surface brightness profile, or that of the underlying bulge 
component. Both have their own merits. The latter pro- 
vides a cleaner measure of the slope of the spheroidal 
component of the bulge, unaffected by inner disks/bars 
or point sources. To date, common practice has been 
to avoid such nuclear features when deriving the inner 
profile slope. On the other hand, a direct measure of 
the slope has the advantage of being independent of any 
model decomposition. Moreover, if the central super- 
massive black holes, known to reside in ellipticals and 
bulges (Kormendy & Richstone 1995) have grown adi- 
abatically within the nuclear star clusters (e.g., Young 
1980; van der Marel 1999), then the profile slope of the 
central cluster is of interest and should not be avoided 
— although resolution does become an issue. 

Several methods have been employed in the past to 
measure the inner profile slope denoted by 7: a single 
power-law (e.g., Phillips et al. 1996); a double power-law 
(Ferrarese et al. 1994); a Nuker fit (L95); a measure of the 
logarithmic derivative of the best-fit Nuker-law at either 
some fixed radius (usually denoted as 7/, e.g., R01, and 
here denoted 7') or over some small interval, e.g. 0.1- 
0.5 arcseconds (usually denoted as (7}; e.g. Stiavelli et 
al. 2001); and the logarithmic derivative of the best-fit 
Sersic model (denoted here as 7'; Graham & Guzman 
2003). For the Sersic model, 7' can easily be evaluated 
at any radius R, such that 

7s '(i?) = -d\ogI{R)/d\ogR = Ji -j , (15) 

where b n w 1.9992n — 0.3271, and R c is the effective 
radius. 

For our sample, we have computed the inner profile 
slope following all of the methods outlined above, except 
for (7), for lack of a precise definition of this quantity. 
The various determinations of the inner profile slopes 
are given in Table O These include direct power-law 
least-squares fits to the profiles over ranges of 0.1 < 
i?/arcsec < 1 and 20 < R/pc < 200; 7^(0.3"); 7 S '(0.3"); 
and, 7 S '(0.15 R c ). Comparison of the results indicates 
that direct power-law fits over ranges of 0.1 < R < 1 
arcsec or 20 < R/pc < 200 yield slopes 7 quite simi- 
lar to 7^(0.3") derived from a Nuker-fit. Our Nuker fits 
extended inward to R = 0.1", hence the similarity of 
7/(0.3") to a direct power-law fit is not surprising: these 
three fits encompass any compact nuclear components. 
The values of 7 S '(0.3") and 7 s '(0.15i? e ), which come from 
the Sersic fits to the host bulge (§E]), are equal or lower 
than the other values. They provide the slope corre- 
sponding to the underlying bulge components, and, as 
such, they are the most directly comparable to the (7) 
values presented by e.g., Carollo & Stiavelli for late- type 
spirals, and those of Stiavelli (2001) for dwarf ellipticals. 

Figure[jj shows the B-band absolute magnitudes of the 



bulges plotted against the values of 7^ derived at 0.3", 
and of 7' derived at R/R c = 0.15. The B-band absolute 
magnitudes of the bulges were derived from the galaxy 
total corrected B-band magnitudes from the RC3 and 
the bulge-to-disk ratios derived in this paper; this pro- 
cedure avoids £?-band bulge-disk decompositions which 
are heavily affected by dust. Also included in the figure 
are the elliptical galaxy data from L95 (Nuker-model 7) 
and the dwarf elliptical data from Stiavelli et al. (2001; 
mean slope (7} between 0.1" and 0.5") and Graham & 
Guzman (2003; 7' at R = 0.2" from Sersic fit). Inner 
slopes for bulges cluster above 7 = 0.5 when the light 
of nuclear components is included (top panel of Fig. [7]), 
and matches that of ellipticals of similar luminosities. 
When nuclear components are excluded (bottom panel 
of Fig. 0), slopes trace a continuous distribution in the 
range < 7 < 1. 

In a previous study of inner profile slopes of galaxy 
bulges, Carollo & Stiavelli (1998) obtain a similar range 
of < 7 < 1 values to us. They however differ in the in- 
terpretation of the results. Because their distribution of 
(7} shows two distinct clumps, they argue that two sepa- 
rate families of bulges exist, i.e. "exponential" ((7) ^0.3) 
and "14" ((7} ~ 0.8) bulges. Figure[7J shows that our 
bulges do not cluster in two clumps, instead they cover 
a continuous range from the low to the high values of 
the Carollo & Stiavelli distribution. The origin for their 
bi-modal distribution of (7} is unclear. Their (7) val- 
ues come from Nuker fits and not from exponential or 14 
fits, hence the bi- modality should not be a consequence 
of their splitting the sample into " exponential" and " 14" 
classes. But two aspects of their analysis could lead to 
a polarization of their (7) toward high and low values. 
One of these is the subjective choice of the radial range 
for the Nuker fits. While nuclear components are easy 
to identify by eye when the underlying profile is shal- 
low, such as in an exponential bulge, for galaxies with 
steeper underlying bulge profiles, (n > 2), the 'break' 
from the bulge to the nuclear component becomes weaker 
and easier to miss by eye, hence the radial range occupied 
by the nuclear component can easily be included when 
performing the Nuker fit, yielding a biased, higher value 
of (7) than that corresponding to the underlying bulge. 
Simultaneous fits to the bulge and the nuclear compo- 
nent are needed in those cases to remove such bias from 
the determination of the bulge profile slope. Another 
aspect that may polarize the (7} distribution into high 
and low clumps may be the derivation of bulge absolute 
magnitudes without a bulge-disk decomposition, and ap- 
pying exponential- or 14-constrained fits to profiles from 
small-field HST/WFPC2 images: 14 models are known 
to over-estimate the flux of R x / n systems if n < 4, while, 
exponential models are known to underestimate the flux 
in R 1 / 11 systems which have n > 1 (e.g., Graham 2001a). 
Finally, sample selection may also be important. Their 
sample includes many barred galaxies, and is overall later 
type than ours, though we note that many of our bulges 
show exponential-like profiles, i.e., profiles which Car- 
ollo & Stiavelli associate with the low- (7) clump. The 
above arguments, together with the distributions shown 
in Figure suggest that the slopes of the total galaxy 
profiles (including nuclear components), cluster around 
0.5 < 7 < 1, but bulge profiles as a class cover a continu- 
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ous distribution of nuclear slopes in the range < 7 < 1. 

6. CONCLUSIONS 

At HST resolution, nuclear photometric components, 
in addition to the Sersic bulge and the exponential 
outer disk, are exceedingly common (~90%) in early- to 
intermediate-type disk galaxies. Spatially-resolved nu- 
clear components are found in 58% of our sample. These 
components are geometrically flat systems, and could be 
disks, bars or rings. The ones detected have comparable 
central surface brightness to the underlying bulges, but 
fainter such systems may exist. The isophotal signatures 
indicate total sizes of a few hundred pc, similar to those 
of inner bars in double-barred galaxies. Often, such com- 
ponents are reddened by dust; the evidence from optical 
and NIR colors, presented in Paper I, as well as their 
high densities, suggest that they are old rather than late 
additions to the bulges. 

A majority of the galaxies (~58%) harbor sources 
unresolved by #5T/NICMOS2. They are most 
likely star clusters, with luminosities corresponding 
to 10-20 globular clusters, and similar to other 
unresolved sources found in the nuclei of ellipti- 
cals, dwarf ellipticals and bulges. When com- 
bined with similar nuclear components in dE galax- 
ies, their photometric masses scale with spheroid mass 
as Mp S /M q = 10 7 - 73±016 (Xbui/10 10 M q )°- 76±013 . 
Our central star clusters fall above the faint-ward ex- 
trapolation of A'l.-A^Bui relations derived by HR04 
or F06. In order to extend a CMO-style relation to 



faint spheroid luminosities, a moderate non-linearity is 
needed, and we propose the relation A4cmo/A^0 = 
1O 7 - 51±o - O6 (7Wbui/1O 1o X )°- 84±0 06 . But we see addi- 
tional difficulties with the CMO picture in that all of our 
PS show masses above the cluster-black hole transitional 
mass of 10 7 A / f© proposed by Wehner & Harris (2006). 

Bulge surface brightness profiles rise inward to the 
limit of the #ST/NICMOS resolution, -10 pc for the 
current sample. While the inner bulge profiles devi- 
ate from pure power-laws, "break radii" in a Nuker-law 
sense are not present. Structurally, the bulges of early- 
to intermediate-type galaxies may be globally grouped 
with the " power-law" , intermediate- and low- luminosity 
elliptical galaxies. Negative logarithmic nuclear profile 
slopes of the Sersic bulge components 7 cover a continu- 
ous range of < 7 < 1, overlapping with dwarf-ellipticals 
at the faint end, and with intermediate luminosity ellip- 
ticals at the bright end. We find no evidence to support 
a bimodal distribution of 7 reported by others. 
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Space Administration. This research has made use of 
the HyperLeda database. The United Kingdom Infrared 
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behalf of the U.K. Particle Physics and Astronomy Re- 
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APPENDIX 

SURFACE BRIGHTNESS, ISOPHOTAL, DYNAMICAL, AND COLOR PROFILES 

Figure KTl shows the surface brightness profiles together with dynamical, isophotal and color profiles for each galaxy. 
The top panel shows the combined HST+GB profile, and Sersic+exponcntial fit carried out excluding the inner 
1 arcsecond. Beneath the galaxy names is the fit type code corresponding to the full fit including additional nuclear 
components, from Table [2] The second panel gives the residuals from the fit. The third panel shows the minor-axis 
velocity dispersion profiles, folded around the origin, from Falcon-Barroso et al. (2003); the abscissa has been scaled 
to the geometric mean radius using the mean of the ellipticity profile. The fourth and fifth panels show the combined 
HST+GB minor-axis B — I and I — H color profiles on the side of the bulge not seen through the disk (the "dust-free" 
side), from Paper I; the abscissa has been scaled to the geometric mean radius using the mean of the ellipticity profile. 
The sixth panel shows the ellipticity profile, and the bottom panel shows the 4th order Fourier cosine term c4, from 
the ellipse fits; inside 6 arcsec {triangles) the latter two profiles are derived from HST/NICMOS F160W images, while 
outside 6 arcsec (squares) they are derived from the UKIRT if-band images (Peletier & Balcells 1997). 
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Fig. Al. — Profiles of surface brightness, velocity dispersion, B — I , I — H , ellipticity, and c4 diskincss/boxincss coefficient. Abscissa 
is the logarithm of the geometric mean axis in arcsec. fin- combined HST+GB surface brightness profile, and model fits excluding the 
central arcsec (dashed line: Sersic bulge; dotted line: exponential disk). A/j,h'- residuals from the fit (data minus model). <r: minor-axis 
velocity dispersion profiles from Falcon- Barroso et al. (2003), folded around the origin (triangles: dust-free side; inverted triangles: side 
seen through the disk; abscissa has been projected to the geometric mean radius using the mean of the ellipticity profile; the central velocity 
dispersion measurement has been plotted at r = 0.18 arcsec, corresponding to one half pixel of the spectrograph.) B — I, I — H: minor-axis 
color profiles for the dust-free side of the bulge, derived from combined HST+GB profiles, from Paper I. Abscissa has been projected to 
the geometric mean radius using the ellipticity profiles, e: triangles arc the ellipticity profile from ellipse fits to the HST/NICMOS F160W 
images; squares are the e profile from ellipse fits to GB iiT-band images, from Peletier & Balcells (1997); GB e values are generally kept 
fixed in the outer, low-S/N region of the images. c4: fourth-order cosine term of the residuals from ellipse fits (triangles: HST F160W 
images; squares: GB iiT-band images). 
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Fig. Al. — continued. 
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